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Detecting chiral singularities in lattice QCD at strong coupling 

Costas G. Strouthos a 

a Division of Science and Engineering, Frederick Institute of Technology, Nicosia 1303, Cyprus 

We study the difficulties associated with detecting chiral singularities in strongly coupled lattice QCD at 
fixed nonzero tempearture. We show that the behavior of the chiral condensate, the pion mass and the pion 
decay constant, for small masses, are all consistent with the predictions of Chiral Perturbation Theory (ChPT). 
However, the values of the quark masses that we need to demonstrate this are much smaller than those being 
used in dynamical QCD simulations. 
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1. Introduction 

Given the difficulties associated with under- 
standing chiral singularities in a realistic QCD 
calculation, we explore this subject in strong cou- 
pling lattice QCD with one staggered fermion at 
finite temperature. The model undergoes a sec- 
ond order transition which belongs to the 0(2) 
universality class [Q. At fixed T < T c , within the 
0(2) scaling window, the long distance physics 
of our model is described by a 3d 0(2) contin- 
uum scalar field theory in its broken phase. At 
this temperature there is a range of quark masses 
where the light pions are describable by a 3d con- 
tinuum chiral perturbation theory (ChPT). The 
effective action can be written in terms of 5, an 
0(2) vector field with the constraint S ■ S = 1. 
At the lowest order this is given by 



d A x 



F 2 



U) 



where h is the magnetic field. F is the pion decay 
constant and S is the chiral condensate. In three 
dimensions F 2 has dimensions of mass. 

We connect our model to ChPT and minimize 
lattice artifacts by choosing T and the quark mass 
m such that 1 <C £ <C 1/M„, where M n is the pion 
mass. In this region we expect that S, M T and 
F 2 satisfy the expansion 



(O) = zq + Z\\fm + z 2 m + z 3 m\/m + .. 



(2) 



where the ^pm behavior is the power-like singu- 
larity arising due to the infrared pion physics P] ■ 



The partition function of lattice QCD with one 
staggered fermion interacting with U(3) gauge 
fields can be mapped into a monomer-dimer sys- 
tem in the strong coupling limit 3 . In this study 
we use a very efficient algorithm discovered re- 
cently to solve the model in the chiral limit 0]. 
In this work we choose L x = L y = L. One 
can study the thermodynamics of the model by 
working on asymmetric and anisotropic lattices 
with L t <C L and allowing the temporal staggered 
fermion phase factor to vary continuously p]. 

In order to study chiral physics we focus on the 
chiral condensate S m , the chiral susceptibility %, 
the helicity modulus Y and the pion mass M„ . As 
discussed in jS] , one can define the pion decay con- 
stant at a quark mass m to be F^ = liniL^oo Y. 
For m = we then obtain F = lim m ^o F m , the 
pion decay constant introduced in Eq.QJ. We 
can also define £ = lim m _>o E m , where E m is the 
infinite volume chiral condensate. 



2. Results 

For massless quarks we fit F 2 vs. T to F 2 = 
C(T C ~ T) v for T < T c . In the fit we fix T c = 
7.47739(3) [TJ and v = 0.67155(27) 0. For 7.05 < 
T < 7.42 we get C = 0.2217(1) with X 2 /DOF = 
1.03. Including T = 7.0 in the fit makes \ 2 /DOF 
jump to 2.0. At T close to T c , £ is very large and it 
is difficult to satisfy -jg- ^> £. On the other hand 
when £ ~ 1 lattice artifacts become important. 
We estimate that T = 7.0 is at the edge of the 
0(2) scaling window. All the results discussed 
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Figure 1. Plot of F m vs. m. The solid line rep- 
resents the function F m = 0.36747[1 + 34.3m — 
185m^/mJ. 



below were obtained at T = 7.0. We fix L t = 4 
and vary spatial volume from L = 8 to L = 96 
for < m < 0.025. 

The first four terms in the chiral expansion of 
F m in 3d O(N) ChPT are given by ©: 



F m = F[l + ai^/m, + a 2 m + a^m\fm\. 



(3) 



where a\ cx (N — 2). Since in our case N = 2 
we expect ai = 0. We fit data in the range < 
m < 0.00625 (see FigHJ). We find F = 0.3674(1), 
ai = 0.01(2), a 2 = 34.0(7) and a 3 = -182(6) 
withx 2 /DO F = 1-2. The prediction of ChPT that 
ai = is in excellent agreement with our results. 

ChPT predicts that the first four terms in the 
chiral expansion of S m are given by 



S[l + hy/m + b 2 m + b 3 my/m\. 



(4) 



As shown in |B] 6i oc (N—l), implying that b\ ^ 
in our case. At m = we compute S by using 
the finite size scaling formula for \ given by [H] 
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1 + 0.226(7V-1)^ + - 



a 



-(5) 



where a depends on the higher order low energy 
constants of the chiral Lagrangian. When F is 
fixed to 0.3675 obtained earlier, our results for 
16 < L < 144 fit very well to this formula. We 
find S = 2.2648(10), a = 4.6(3) with x 2 /DOF= 



Figure 2. Plot of S m vs. m. The solid line rep- 
resents the function S m = 2.2642[1 + 1.36^/771-1- 
23m—135my / m\. The dashed line in the inset rep- 
resents E m - 2.2642 = 2.2643[1.44vm + 15.7m]. 
For clarity some of the data points are not shown 
in the inset. 



0.87. Fitting our data to Eq.QJ in the region < 
m < 0.00625 gives E = 2.2642(10), b x = 1.36(4), 
b 2 = 23(1), 6 3 = -135(10) with a x 2 /DOF= 1.1. 
We note that we cannot find a mass range within 
our results in which we can find a good fit when 
we fix 62 = &3 = 0. However, in the range < 
m < 0.00175 we can set 63 = to obtain £ = 
2.2643(10), 61 = 1.44(3) and b 2 = 15.2(7) with 
a x 2 /DOF= 1. Note that b 2 changes by about 
30% when this different fitting procedure is used, 
while 61 is more stable. The data and the fits are 
shown in Fig [3 

The first four terms in the chiral behavior of 
the pion mass are predicted to be of the form 



1 + civm + c 2 m + C3my/m 



.(6) 



Further, the chiral Ward identities imply that 
ci = &i/2 — ai. Fitting our data to Eq.© af- 
ter fixing S = 2.2648 and F = 0.3675 obtained 
above, we find c\ = 0.55(16), c 2 = —21(6) and 
c 3 = 110(50) with x 2 /DOF= 0.5. We see that our 
data is consistent with the relation c\ = b\/2 — ai 
although the error in c\ is large. Fixing c\ = 
bi/2 — 0.68 obtained from fitting the chiral con- 
densate, yields c 2 = -25.8(8) and c 3 = 149(14) 
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Figure 3. Plot of M„ vs. m. The solid line repre- 
sents the function M n = 4.095y / m[l + 0.68^/m — 
25.8m + U9my/m\. 



without changing the quality of the fit. This fit 
along with our data for M„ are shown in Fig|3| 

3. Discussion 

It is expected that the natural expansion pa- 
rameter for ChPT in three dimensions is x = 
I '(4ttF 2 ). Let us find the values of x where 
the chiral expansion up to a certain power of x 
is sufficient to describe the data reasonably. The 
solid lines in Fig0| are the best fits to the form 
zq + z\x + zix 1 + Z3X 3 , while the dot-dashed lines 
are the best fits to a smaller range in x with Z3 set 
to 0. As the graph indicates the relative error in 
S m due to the absence of the Z3 term is smaller 
from the error in at a given value of x. In 
order to determine within say 5% we need 
the Z3 term even at x ~ 0.15. Thus, we conclude 
that in our model the chiral expansion converges 
rather poorly, especially for the pion decay con- 
stant. It is striking that almost all previous re- 
liable unquenched lattice QCD simulations with 
staggered fermions have used m > 0.005 (in lat- 
tice units), irrespective of whether the study was 
of QCD thermodynamics or the QCD vacuum, 
whether the studies were done at strong couplings 
or weaker couplings. The value m ~ 0.005 cor- 
responds to 1 ~ 0.18 in our model. We can ask 
how the chiral extrapolations will look if we use 



Figure 4. Plot of F^ and S m vs. x = 
M^/i^nF 2 ). The solid lines and dot-dashed lines 
are fits to two different orders of chiral expansion 
discussed in the text. The dashed line represent 
fits using only the data with x > 0.18. 



the data for x > 0.18 to fit to the expected chi- 
ral form. In Fig0]we show the best fits for large 
masses using the dashed lines. As expected in or- 
der to use chiral extrapolations reliably it is im- 
portant to know the range of their validity. 
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